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ABSTRACT 

We show that observations of pulsars with pulsed optical emission indicate 
that the peak flux scales according to the magnetic field strength at the light 
cylinder. The derived relationships indicate that the emission mechanism is com- 
mon across all of the observed pulsars with periods ranging from 33ms to 385 ms 
and ages of 1000-300,000 years. It is noted that similar trends exist for 7 ray pul- 
sars. Furthermore the model proposed by Pacini (1971) and developed by Pacini 
and Salvati (1983,1987) still has validity and gives an adequate explanation of 
the optical phenomena. 

1. Introduction 

Since the first optical observations of the Crab pulsar in the late 1960s ((4)), 7 more 
pulsars have been observed in the optical bandpass. Of these only 4 have been seen to 
pulsate optically (PSR0540-69, (20); Vela, (37); PSR 0656+14, (32); Geminga, (33)). Whilst 
the three remaining pulsars (PSR0950+08, PSR 1929+10 and PSR1055-52) are thought to 
be mainly thermal emitters in the optical regime, the pulsating objects are either purely 
magnetospheric emitters or a mixture of thermal and non-thermal sources. Two of these five 
pulsars (PSR0540-69 and the Crab) are certainly too distant to have any detectable optical 
thermal emission using currently available technologies. For those with distances less than 
a kiloparsec (Geminga, Vela and PSR0656+14) only the Vela pulsar is too bright for its 
emission to be anything other than non-thermal in origin. For the remaining two pulsars 
(PSR0656+14 and Geminga) their optical emission is thought to be a mix of thermal and 
non-thermal fluxes ((8) and (18)). However for the faintest of the two, Geminga, there is 
disagreement about the level of thermal emission. Broad-band photometric studies have 
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been interpreted as being consistent with the Rayleigh- Jeans tail combined with a cyclotron 
emission feature ((21)). These authors suggest that the optical emission is predominantly 
thermal with an embedded ion cyclotron resonance feature at about 5500 A. More recently, 
(14), derive a complex phenomenological model for this feature coming from a thin plasma 
above the pulsar's polar cap. Their derived surface magnetic field using this model is 
consistent with the independently determined 'spin-down' field {B oc (PP)°-^). Conversely, 
spectroscopic studies have indicated that the emission is a mix of thermal and non-thermal 
radiation with the non-thermal dominating in the optical region ((18)). These authors do 
not see the proposed emission feature, but, there is a statistically weak absorption feature 
at 6400A. Temporal studies of the pulsars indicate that for the Crab and Vela pulsars the 
peak of 7-ray and optical pulses are coincident. A similar coincidence is noted for the 
Geminga pulsar although the optical signal is of lower significance. This can interpreted 
as indicating a common source location for both 7 and optical photons. Recent studies of 
the plateau of the Crab main pulse in the optical indicate that it is of limited (< 15 kms) 
extent ((9)). 

Many suggestions have been made concerning the non-thermal optical emission process 
for these young and middle-aged pulsars. Despite many years of detailed theoretical studies 
and more recently limited numerical simulations, no convincing models have been derived 
which explain all of the high energy properties. There are similar problems in the radio 
but as the emission mechanism is radically different (being coherent) only the high energy 
emission will be considered here. 

In recent years a number of groups have carried out detailed simulations of the various 
high-energy emission processes. These models divide into two broad groups - those with 
acceleration and emission low in the magnetosphere (polar cap models) and those with the 
acceleration nearer to the light cylinder (outer-gap models). Both models have problems 
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explaining the observed features of the hmited selection of high energy emitters. Both 
models also suffer from arbitrary assumptions in terms of the sustainability of the outer-gap 
and the orientation of the pulsar's magnetic field to both the observers line of sight and the 
rotation axis. Furthermore some observational evidence, see for example (6), severely limits 
the applicability of the outer-gap to the emission from the Crab pulsar. However these 
models have their successes - the total polar-cap emission can be understood in terms of 
the Goldreich- Julian current ((12)) from in or around the cap; the Crab polarisation sweep 
is accurately reproduced by an outer-gap variant of (30). However, the most successful 
model which adequately explains most of the high-energy phenomena (both in terms of 
its elegance and longevity) has been proposed by (26) and in modified form by Pacini 
and Salvati (1983 PS83 & 1987 PS87 hereafter). In general they proposed that the high 
energy emission comes from relativistic electrons radiating via synchrotron processes in 
the outer regions of the magnetosphere and that the luminosity is a strong function of the 
period (oc P^^o in the original formulation). In this paper we examine the validity of their 
approach and show that it still adequately explains the observed phenomena. 

It is the failure of the detailed models to explain the high energy emission that has 
prompted this work. We have taken a phenomenological approach to test whether Pacini 
type scaling is still applicable. Our approach has been to try to restrict the effects of 
geometry by taking the peak luminosity as a scaling parameter rather than the total 
luminosity. In this regard we are removing the duty cycle term from PS87. It is our 
opinion that to first order the peak emission represents the local power density along 
the observer's line of site and hence reflects more accurately emission processes within a 
pulsar's magnetosphere. Previous work in this area (see for example (11)) looked at the 
total efficiency, spectral index and found no reasonable correlation with standard pulsar 
parameters - age, period and spin down rate. Their work was hampered by not including 
geometry and being restricted to the then three known pulsed optical emitters. Since then 
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the number of pulsars with observed pulsed optical emission has increased to five. 

2. The Phenomenology of Magnetospheric Emission 

The three optically brightest pulsars (Crab, Vela and PSR 0540-69) are also amongst 
the youngest. However all these pulsars have very different pulse shapes resulting in a very 
different ratio between the integrated flux and the peak flux. In this work we will use the 
peak emission as the primary flux parameter. A number of deflnitions can exist for this and 
in this context we have taken the 95%-95% level of the primary pulse. To flrst order, this 
correlates well with the luminosity per pulse divided by the Full- Width at Half-Maximum 
(FWHM). For the Crab pulsar the cusp is effectively flat over this region ((10)). The 
FWHM can be considered to scale with the pulsar duty cycle. Our proposition is that the 
peak flux represents the local power density within the emitting region with minimal effects 
from geometrical considerations such as observer line of sight and magnetic and rotational 
axis orientation. Table 1 shows the basic parameters for these objects including the peak 
emission (taken as the emission in the 95% - 95% portion of the largest peak). Their 
distances imply that the thermal emission should be low, in all cases < 1% of the observed 
emission, and any such contribution is negligible. 

Of all the optical pulsars Geminga is perhaps the most controversial. Early observations 
((16), and (2)) indicated that Geminga was an ^ 25.5 my object. Subsequent observations 
including HST photometry appeared to support a thermal origin for the optical photons, 
albeit requiring a cyclotron resonance feature in the optical ((21)). The high-speed optical 
observations of (33) combined with spectroscopic observations ((18)) contrast with this 
view. Figures 1 and 2 show how this apparent contradiction could have arisen. Figure 1, 
based upon data from (21) shows the integrated photometry. It is possible to flt a black 
body Rayleigh- Jeans tail through this, but only with the fltting of a cyclotron resonance 
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emission feature at about 5500 A. Figure 2 however shows the same points plotted on top 
of the Martin et al spectra, where we have also included the pulsed B point ((33)). All the 
data points are consistent within their error bars. This combined data set indicates a fairly 
steep spectrum (with spectral index of 1.9) consistent with magnetospheric emission and a 
weak thermal component, without the requirement for a cyclotron feature. It was on the 
basis of these results that (8) were able to give an upper limit of Rqo of about 10km, by 
considering the upper limits to the unpulsed fraction of the optical emission from Geminga 
as an upper limit for the thermal emission. A simpler view is that the soft X-ray and 
EUV data is predominantly thermal emission with a magnetospheric component becoming 
dominant at about 3500 A. Figure 3 illustrates this with the thermal and magnetospheric 
components from (8) combined with the broad-band points of (21) and (33). We should 
also note that the optical pulses seen by (33) are coincident with the observed 7 ray peaks. 
Clearly more spectroscopic data, particularly if temporally resolved, will be crucial to 
determining the exact mix of thermal and non-thermal emission. 

As regards PSR 0656+14 the pulsar is generally agreed to be predominately a non- 
thermal emitter in the optical, becoming predominantly a thermal emitter at wavelengths 
shorter than about 3000A((25)). However, there is a discrepancy between the radio 
distance based upon the dispersion measure and the best fits to the X-ray data. From 
radio dispersion measure a distance of 760 ± 190pc can be derived, at odds with the X-ray 
distance of 250 — 280pc from A'^^ galactic models. Clearly more observations are needed 
to determine a parallax - both radio and optical. This discrepancy in distance leads to an 
ambiguity in the total luminosity. For this paper we have taken the lower distance measure. 

In Table 1 we have indicated the peak fiux normalised to the Crab pulsar. When 
comparing the individual pulsars, which will have a range of viewing angles and differing 
magnetic/rotation axes, we have to account for the viewing angle (related to the pulse 
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duty cycle and separation) as well as the total flux. Furthermore for each source, we have 
to account for the shape of the pulsar's spectrum. For all the observed pulsars, with 
the exception of Vela, the low energy cut-off is above 7000 A((23)). PS87 indicates for 
emission above the low energy cut-off that the ratio of the fluxes from two different pulsars 
can be given by Equation 1 if one ignores the effect of duty cycle and pitch angle (the 
suffices refer to each pulsar). Here Fi^n refers to the flux at the observed frequency z/ for 
pulsar n. Similarly for the magnetic field strength, B, and period, P. The observed energy 
spectrum exponent is given by a„. The duty cycle can be accounted for by only considering 
the peak emission. The pitch angle is beyond the scope of this work and assumed to 
first order to be invariant. Equation 2 shows the the same formulation for the outer field case. 



Fv,l V Bifl Pi 

Scaling to the transverse field would give:- 

Fu,i V -Bi,o -Pi 

Given the observed peak luminosities we investigated the correlation between the peak 
emission and the surface field and the tangential light cylinder field. Figure 4 shows the 
relationship between the peak luminosity and the outer magnetic field, S-r, Goldreich-Julian 
current and canonical age (^)- A clear correlation is seen with all these parameters. In 
this paper we investigate the implications of a correlation between the peak luminosity and 
the transverse field. We accept that the strong correlation with G-J current would underpin 
both emission from polar cap as well as outer regions. 

A weighted regression of the form :- 



Peak Luminosity oc Bj, 

was performed for the empirical peak luminosity leading to a relationship of the form :- 

Peak Luminosity oc B^^^^'^-^^ 

which is significant at the 99.5% level. Figure 5 shows the predicted peak luminosity 
from Equation 2 against our observed peak luminosity accounting for the differing observed 
energy spectrum exponent at 4500 A. The slope is 0.95 ± 0.04 and significant at the 99 
% level. Importantly, we note that the fiattening of the peak luminosity relationship with 
the outer field strength for the older, slower pulsars, see Figure 4, is consistent with their 
having a significantly steeper energy spectrum than the younger pulsars, see Table 1. 
Whilst, from PS87, it is possible that this indicates that the emission zone is optically thick 
at these frequencies, alternatively it might refiect a larger emission region for these pulsars 
compared to the younger ones. We also note that the generation parameter concept ((38)) 
indicates that the 7-ray spectral index scales with the average number of generations within 
the e~^ /e^ cascade. Our observed correlation between peak luminosity and age (which 
linearly scales with the generation parameter) and similarly between spectral index and age 
indicates a link between optical luminosity and 7-ray luminosity phenomena. 

We note as well as that similar trends can be seen in 7-rays. Figure 6 shows the 
correlation between the peak 7 emission as a function of transverse field. This indicates a 
regression of the form 

Peak 7 Luminosity oc 5°-®^^°-^^ 

significant at the 99.3 % level, consistent with the observed steeper distribution 
seen in high energy 7-rays - see Table 2 and Equation 2. 

It seems clear from both polarisation studies ((34); (30)) and from this work that we 



-9- 

expect that the optical emission zone is sited towards the outer magnetosphere. Timing 
studies of the size of the Crab pulse plateau indicates a restricted emission volume (~ 15 
kms in lateral extent) ((9)). Importantly, the simple relationships we have derived here 
indicate that there is no need to invoke complex models for this high energy emission. 
Observed variations in spectral index, pulse shape and polarisation can be understood 
in terms of geometrical and absorption factors rather than differences in the production 
mechanism. 

Combining our results and those of (11), we can begin to understand how the high 
energy emission process ages. Goldoni et al compared the known spectral indices and 
efficiencies in both the optical and 7-ray regions. They noted that the spectral index 
flattened with age for the 7-ray pulsars whilst the reverse was true for the optically emitting 
systems. They also noted a similar trend reversal for the efficiency, with the 7-ray pulsars 
becoming more efficient with age. We note (see the bottom panel in Figures 4 and 6) a 
similar behaviour with the peak emission. From the temporal coincidence between the 7-ray 
and optical pulses, it seems likely that the source location is similar for both mechanisms. 
One explanation is that we have a position where by from the same electron population 
there are two emission processes - expected if we have curvature for the 7-ray photons 
and synchrotron for the optical ones. It seems likely that the optical photon spectrum has 
been further modified to produce the reversal in spectral index with age. The region over 
which the scattering can occur would scale with the size of the magnetosphere and hence 
with age. With the outer magnetosphere fields for these pulsars being < 10^ G, electron 
cyclotron scattering is not an option. However synchrotron self- absorbt ion could explain 
the observed features. In essence we would expect the most marked flattening to be for the 
Crab pulsar, where the outer fleld strengths are of order 10^ G, and less so for the slower 
and older systems. 
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These results (both optical and 7-rays) are consistent with a model where the 7 and 
optical emission is coming from the last open-field line at some constant fraction of the 
light cylinder. The drop in efficiency with age for the production of optical photons points 
towards an absorbing process in the outer magnetosphere. Clearly more optical and 7-ray 
observations are needed to confirm these trends. 

3. Conclusion 

We have shown that the peak optical luminosity is the key important parameter that 
scales with other observed pulsar attributes. We note that similar behaviour can be seen 
in the 7-ray regime. To first order we confirm that the model first proposed by Pacini 
in 1971 still has validity. There is a proviso that there is a strong correlation with the 
Goldreich- Julian current which underpins both Pacini scaling and other models. From a 
more detailed anaysis of the EGRET database (19) showed a similar functional form to 
our derived relationship for the total luminosity against period and surface magnetic field. 
We also show in Table 3 what the expected luminosity would be from a number of X and 
7-ray emitting pulsars. We have chosen these on the basis of their known duty cycles in the 
high energy regime. We have also considered the expected fiux from the anomalous X-ray 
pulsars, a class of magnetars. If the observed P and P relationship can be interpreted in 
same way as for normal pulsars to derive a canonical surface field then we can estimate 
what the expected luminosity of these object would be. The derived luminosity is very low 
- a reflection of the weak light cylinder field. A recent VLT observation ((15)) has found no 
optical candidates down to 25.5 R magnitude - consistent with no magnetospheric emission 
and predominantly thermal output. Problems remain however as to the source of the 
thermal emission. We note that, in the outer magnetosphere of these objects, the magnetic 
field is far too weak for optical non-thermal emission processes. 
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Of crucial importance in the future will be the determination of the low-energy cutoff 
and the polarisation sweep through the optical pulse. Also of interest will be the shape 
of the pulse - in particular the size of any plateau which scales as the size of the emission 
zone. Optical observations of millisecond pulsars will also be important, to see if they 
follow the age or field trends noted here. Of interest in this regard is the result from X-ray 
observations which indicate that the luminosity of millisecond pulsars scales in a similar 
fashion to normal pulsars ((1)). All of these parameters are measurable with existing and 
emerging technologies and telescopes. 

Finally, the recent detection of a 16 ms pulsar in the LMC (PSR J0537-6910 ((17)) 
which has defied optical detection despite its low period can possibly be understood in 
terms of its age (~ 5,000 years), which on the basis of Figure 3 indicates peak luminosity a 
factor of ~ 10^ times lower than the Crab. This is at the limit of some of the recent optical 
searches - (22) and (13). 
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Fig. 1. — Photometry of Geminga. The integrated photometry and the thermal fit to the 
ROSAT X-Ray and EUVE data are both shown. 
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Fig. 2. — Spectra and Photometry of Geminga. The spectrum taken from (18) (sohd hne) and 
1 a error hmits (dotted hne) are shown. Note the agreement between it and the integrated 
photometry from Figure 1. B^^^^ is the pulsed flux from (33) and other points are from 
(21) and (3). 
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Fig. 3. — Derived thermal and non-thermal components for Geminga. Also plotted are the 
broad-band points points from figures 1 and 2. 
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Fig. 4. — Peak optical luminosity vs light cylinder field, Goldreich & Julian current and 
canonical age. Also shown is the efficiency of the peak emission against age. The peak 
luminosity has been normalised to the Crab pulsar. The error bars represent both statistical 
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Fig. 5. — Predicted peak optical luminosity from equation 2 versus observed peak emission. 
Also shown is the weighted fit described in the text. The weighting was based upon the 
observational uncertainities in flux and distance. 
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Fig. 6. — Peak 7-ray luminosity vs light cylinder field, Goldreich & Julian current and 
canonical age. Also shown is the efficiency of the peak emission against age. The peak 
luminosity has been normalised to the Crab pulsar. The 7-ray peak luminosity has been 
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Table 1: Main characteristics of Optical Pulsars: B5 & B^c", the canonical surface and 
transverse magnetic field at the light cylinder respectively; Int and Peak refer to the optical 
luminosity, integrated and peak, at the indicated distance in the B band. 
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2 


33 


42 


3.09 


12.6 


6.1 


10« 


10" 




-0.11 




15000{?) 


Vela 


0.5 


89 


11 


4.11 


12.5 


4.8 


27 


21 




0.2 




6500(?) 


PSR0545-69 


49 


50 


40 


3.20 


12.7 


5.7 


1.1 10" 


1.4 lO'"" 




0.2 




>7000 


PSR0656+14 


0.25(?) 


385 


1.2 


5.50 


12.7 


3.0 


1.8 


0.3 




1.3 




>8000 


PSR0633+17 


0.16 


237 


1.2 


4.99 


12.2 


3.2 


0.3 


0.1 




1.9 




>8000 



-22- 



Table 2: Main characteristics of 7-ray Pulsars: B^ & B^c the canonical surface and transverse 
magnetic field at the light cylinder respectively; 7 Int and Peak refer to the integrated and 
peak 7-ray luminosity at the indicated distance for E > 100 MeV. 



Name 




D 


P 




P 




Bs 


Blc 


Int Lumin. 


Peak Lumin. 


Spectral Index 






(kpc) 


(ms) 


10 


-14 ^ 


,/s 


log(G) 


log (G) 


Crab = l 


Crab=l 






Crab 




2 


33 




42.1 




12.58 


5.97 


1 


1 




2.15 


Vela 




0.5 


89 




12.5 




12.53 


4.64 


0.0480 


0.0299 




1.70 


PSRK 


)55-52 


0.5-1.5 


197 




0.6 




12.03 


3.11 


0.0124 


0.0037 




1.18 


PSRl- 


r06-44 


2.4 


102 




9.3 




12.49 


4.42 


0.1380 


0.0368 




1.72 


PSROe 


533+17 


0.16 


237 




1.1 




12.21 


3.05 


0.0019 


0.0008 




1.50 


PSRl! 


)51+32 


2.5 


40 




0.6 




11.69 


4.86 


0.0500 


0.0187 




1.74 



-23- 



Table 3: Predicted Luminosity of X and 7-ray emitting pulsars. The duty cycle has been 
estimated from 7-ray observations. Also included are the nearby millisecond pulsar PSR 
J2322+2057 ((24)) and the anomolous X-ray pulsar 1E1841-045 ((Vasisht & Gotthelf)). 



Name 


D 


P 


P 


Bs 


Blc 


Duty 


Prcd. Lumin. 




(kpc) 


(ms) 


10-1* J./J. 


log(G) 


log (G) 


Cycle 


fiCrab 


PSR1055-52 


1.5 


197 


0.6 


12.03 


3.11 


0.2 


0.01 


PSR1706-44 


1.8 


102 


9.3 


12.49 


4.42 


0.14 


35 


PSR1951+32 


2.5 


40 


0.6 


11.69 


4.86 


0.08 


670 


PSR1821-24(M28) 


5.1 


3 


1.1-* 


9.3 


5.8 


0.1 


0.3-1 10^ 


PSRJ2322+2057 


0.78 


4.8 


7.0-^ 


8.3 


4.2 


0.3(?) 


8 


1E1841-045 


7 


11770 


4700 


KS 15 


-2 


0.5? 
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